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NONLINEAR ANALYSIS OF MULTILAYERED, 
HORIZONTALLY NAIL-LAMINATED WOOD BEAMS 


by 


Dave R. Bohnhoff, Assistant Professor 
Agricultural Engineering Department 
University of Wisconsin - Madison 
Madison, Wisconsin 


ABSTRACT 


A simple means to analyze multilayered, horizontally, 
nail-laminated beams using the direct stiffness method of 
finite element analysis is presented. Special element 
stiffness matrices and element load modifying vectors are 
formulated for use in modeling nail fasteners. 


INTRODUCTION 


Beams which are fabricated by joining individual pieces of wood togeth- 
er fall into four main categories: 


Horizontally, glue-laminated assemblies. 
Horizontally, mechanically-laminated assemblies. 
. Vertically, glue-laminated assemblies. 

. Vertically, mechanically-laminated assemblies. 


FUNE 


As indicated by the titles of the four preceding groups, laminated 
assemblies are categorized according to the orientation of the assembly 
with respect to the applied load, and according to the method used to 
join the wood pieces together. Briefly, horizontally laminated 
assemblies are members designed to resist loads applied normal to 
interlayer planes. Conversely, vertically laminated assemblies are 
designed to resist loads applied parallel to the planes of contact 
between layers. Figure 1 illustrates the differences between these two 
assembly types. Glue-laminated or glulam members are’assemblies in 
which wood layers are rigidly bonded together with an adhesive. In the 
analysis of these assemblies, it is generally assumed that there is 
complete composite action, that is, that there is no slip between the 
layers when the assembly is loaded. Mechanically laminated assemblies 
- assemblies in which layers are joined together with some type of 
mechanical fastener (nails, bolts, etc.) - differ from glulam members 
in that there can be considerable slip between the individual layers. 


Almost all glulam assemblies are horizontally laminated. This is 
primarily due to the fact that all curved beams are horizontally 
laminated assemblies. It is also important to note that most straight 
glulam beams are also horizontally laminated. This is primarily due to 
the fact that in horizontally laminated assemblies, higher grade lumber 
can be placed farther from the neutral axis and thus is used more 
effectively. Another advantage of horizontally laminated assemblies is 
that variations in the modulus of elasticity of individual lamina occur 
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FIG. 1.- Vertical Versus Horizontal Lamination 


in a direction parallel to the applied load. Consequently, the shear 
center of the assembly is not located off to one side of the beam, and 
therefore, twisting and lateral displacements will be minimal under 
planar loads. The major disadvantage of horizontally glulam assemblies 
is that their strength is highly dependent on the quality of the splice 
joints in the outside laminae. 


Mechanical fasteners are generally used when facilities. for gluing wood 
together are not available, or when it is not economical or practical 
to use glue. Of the different mechanical fasteners available for 
laminating, nails are far and away the most popular. "Nail-laminated" 
assemblies find considerable use as girders and columns in residential 
as well as commercial and industrial buildings. 


Unlike for glulams, in nail-laminated assemblies vertical lamination is 
perferred over horizontal lamination. This is because a measurable 
amount of interlayer slip must occur before nails transfer any 
appreciable amount of load. This slip is critical in horizontally 
laminated assemblies because it can reduce the amount of composite 
action to the point where individual layers virtually work 
independently to resist applied loads. Interlayer slip in vertically 
laminated assemblies is less significant in terms of its affects on 
overall assembly strength and stiffness. Nails in vertically laminated 
assemblies basically act as the load distributing elements in a load 
sharing system. 


The significant effect that interlayer slip has on the stiffness of 
horizontally laminated assemblies has led a number of researchers to 
develop methods for predicting their behavior. The bulk of this 
research has involved floor systems fabricated with elastomeric 
adhesive and nails, and has resulted in the advancement of specialized 
finite element analysis programs. The first of these programs, FEAFLO 
was developed by researchers at Colorado State University (Thompson et 
al, 1975) and was based on earlier work by Goodman (1969). Other large 
programs include: NONFLO - an extension of FEAFLO which accounts for 
nonlinear nail-joint behavior (Wheat et al, 1983), and FAP - a program 
written by Foschi (1982). A simplified method for computing the 
deflections of wood-joist floor systems was presented by McCutcheon in 
1977 and was based on research by Kuenzi and Wilkinson (1971). 
McCutcheon (1984) extended his original analysis to account for 
variable joist stiffness and two-way action. A simple BASIC program 
for implementing the modified procedure is listed in the 1984 
publication. 





Other methods for predicting the stiffness of horizontally laminated 
assemblies with interlayer slip have been developed independently of 
floor system research. Itani and Hiremath (1980) presented a method of 
analysis for what they referred to as layered diaphragms. The 
procedure was implemented in a program called DAD (Difference Analysis 
of Diaphragms). More recently, Kamiya (1985) formulated a procedure 
for predicting the nonlinear bending behavior of 2 and 3 layered 
nail-laminated beams and McCutcheon (1986) developed a simple procedure 
for computing the linear-elastic composite stiffness of "a wood bending 
member with sheathing attached nonrigidly to one or both edges." 


It is evident from a review of the literature that most of the 
theoretical analyzes for horizontally laminated assemblies are based on 
the same assumptions. Therefore, the final governing equations for the 
theories presented are essentially the same. The analysis methods 
which currently are available for floor systems can be divided into 
those which are embodied in finite element programs and those which are 
not. For the most part, methods which don’t involve finite element or 
finite difference techniques are more limited in terms of their 
applicability since more restrictions are generally placed on the model 
in an attempt to simplify the solution process. For example, the 
analysis may be limited to 2 or 3 layered assemblies, or only certain 
load configurations or support conditions may apply. 


Finite element methods of analysis have increased in popularity due to 
the advancement of highly reliable and affordable desktop computers. 
Today, one would be hard pressed to find a structural engineer who is 
not familiar with the direct stiffness method of structural analysis, 
or one who does not have access to a simple finite element program for 
analyzing plane-frame structures. At the same time, it is fair to say 
that few engineers have ready access to a program for analyzing 
horizontally nail-laminated assemblies. Unfortunately, what few people 
realize is that it doesn’t take much to turn a conventional program for 
analyzing plane-frame structures into a program for analyzing 
horizontally nail-laminated assemblies. 


In the following paper, a procedure for analyzing multi-layered, 
horizontally nail-laminated assemblies using the direct stiffness 
method of analysis will be presented. Direct physical arguments will 
be used to formulate the stiffness matrix for both a linear and 
nonlinear nail element. This type of stiffness matrix formulation is 
perferred over methods based on the principle of minimum potential 
energy because it is easier for those with only an introductory 
background in structural analysis to conceive. 


CONSTRUCTION OF THE MODEL 


The first step in modeling a multi-layered, horizontally, nail-lami- 

` nated assembly is to properly divide the beam into elements. This is 
accomplished by vertically sectioning the beam at any location along 
the beam where there is: (1) a support, (2) a fastener, (3) an applied 
load, or (4) a member end. Additional sectionings should be made at 
any point where displacements are desired. The next step is to place a 
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node where each section line intersects the neutral axis of a wood 
layer. If there are n wood layers and s section lines this results in 
n*s structural nodes. Figure 2(a) shows the location of the section 
lines for a three layer beam and Figure 2(b) shows the same beam after 
placement of the nodes. 


BEAM SECTIONING Load 
O = Section Line Number 





Nail 





(b) 
SCHEMATIC REPRESENTATION Horizontal Lines Connecting 


Nodes Represent Wood Elements 





Vertical Lines Connecting 
X Nodes Represent Nail Elements 


(c) 


FIG. 2.- Three Layer Assembly. (a) Sectioning. (b) Node Placement. 
(c) Schematic Representation 


BASIC MODELING ELEMENTS 


Ohly two basic elements, a plane-frame wood element and a special nail 
element, are required to model horizontal, nail-laminated beams. As 
shown in Figure 2(c), each of these elements is associated with only 
two structural nodes. The nodes define the ends and neutral axis of 
the plane-frame wood elements. Since there are usually s-1 wood 
elements per lamina, each model is normally comprised of n*(s-1) wood 
elements. It is important to note that because section lines were 


placed wherever there was a nail, each plane-frame wood element is 
assured to be void of nails. 


Nail elements are used to model nail joints. One nail joint is defined 
as the connection between two wood layers. A nail in double shear 
would therefore be modeled using two nail elements. A nail element is 
associated with the nodes on the section line which lie immediately 
above and below the nail joint. 


INDEPENDENT STRUCTURAL DEGREES OF FREEDOM 


Since this is a planar analysis, each node is associated with three 
degrees of freedom, a vertical (Y) translation, a horizontal (X) 
translation, and a rotation (about the Z axis). Of the 3*n degrees of 
freedom at a section line, only n+2 are independent. This is due to 
the fact that all nodes which lie on the same section line are forced 
to have the same Y displacement and consequently, the same Z rotation. 


PLANE-FRAME WOOD ELEMENT STIFFNESS MATRIX 


In the direct stiffness method of analysis, a stiffness matrix is 
systematically formulated for each element and added into the 
structural (or global) stiffness matrix. The stiffness matrix for an 
element relates the displacements of ‘the element nodes to the forces 
which must be applied to the nodes to keep the element in a state of 
static equilibrium. In matrix notation this relationship is written 
as: 


(k]{d} = {r} [1] 


where: 
[k] = element stiffness matrix. 


{d} = vector of element node displacements. 
{r} = vector of element nodal loads. 


Ý Yi» qi vjes] 





“is Pi 


X 8, mi Oj, mj 
FIG. 3.- Plane-Frame Element 


Figure 3 gives the nomenclature and positive sign convention for the 
nodal forces and displacements acting on the plane-frame element 
located in layer k between section lines i and j. Since each of the 
plane-frame wood elements is void of nail joints and assumed to behave 
linearly-elastically, the relationship between the nodal forces and 
nodal displacements is given by the conventional plane-frame element 
stiffness matrix. With respect to the sign convention used in Figure 


3, this relationship is given as: 
AL?/I 0 O -AL?/I 0 0 Uik Pik 
0 12 6L 0O -12 6L vi qi 
EI 0 6L 4L? 0 -6L 2L? Oy my 
L3 = [2] 
-AL?/I 0 0 AL? /I 0 0 ujk Pjk 
(0) -12 -6L 0 12 -6L vj qj 
0 6L 2L? 0 -6L 4L? 8; mj 
where 
E = Modulus of elasticity of the wood. 
L = Length of the element or distance between section lines i 
and j. : 
A = Cross-sectional area of the wood element. 
I = Moment of inertia of the wood element about the z axis. 
uj, = Horizontal or X translation of layer k at section line i. 
vi = Vertical or Y translation of all nodes on section line i. 
Əi = Rotation of all nodes on section line i. 
ujk = Horizontal or X translation of layer k at section line j. 
vj = Vertical or Y translation of all nodes on section line j. 
8; = Rotation of all nodes on section line j. 
Pik = Horizontal load applied to the element at section line i. 
qj = Vertical load applied to the element at section line i. 
mj ~ Moment applied to the element at section line i. 
Pjk 7 Horizontal load applied to the element at section line j. 
qj = Vertical load applied to the element at section line jÍ. 
mj = Moment applied to the element at section line j. 


The preceding plane-frame element stiffness matrix has been derived by 
a number of authors (Wang, 1983; Segerlind, 1984; Cook, 1981; Willems 
and Lucas, 1978) and will not be formulated here. 


FORMULATION OF THE NAIL ELEMENT STIFFNESS MATRIX 


Linear Nail-Joint Load-Slip Behavior 


Figure 4(a) shows a typical nail joint. The location of the nail at 
the interface between the top member T and the bottom member B is 
marked by points Cy and Cy. Both points have the same initial X 
coordinates, however Cy belongs to member T and Ch belongs to member B. 
When members T and B slide relative to each other, Cy is displaced from 
Cp an amount A (Figure 4(b)). The quantity A is defined as the slip of 
the nail joint. 


(a) 


(b) 


(c) 
X 


FIG. 4.- Typical Nail Joint. (a) Initial Geometry. (b) Displaced 
Geometry. (c) Corresponding Nail Shear Force F 





As shown in Figure 4(c), the slip A of the nail joint induces a resist- 
ing shear force F in the nail. The relationship between the force F 
and the slip A is referred to as the nail joint load-slip relationship. 
This relationship is dependent on the physical and mechanical 
properties of the nail and the wood as well as the geometry of the 
joint. It is typically determined by conducting a number of nail joint 
load-slip tests. A typical nail joint load-slip curve is shown in 
Figure 5. As the curve indicates, the load-slip relationship of a nail 
joint is usually highly nonlinear. However, below a certain load limit 
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(e.g. the fastener design load limit), the load-slip relationship for 
the joint is often assumed linear and the shear force F is given as: 


F = kg*A [3] 


where kg is the nail-joint stiffness in the linear-elastic range. 


F 
i tan '(k,) 
A 


Nail Slip 
FIG. 5.- Load-Slip Curve for a Typical Nail Joint 


Nail Shear Force 





(c) (d) 


FIG. 6.- (a) Nail Shear Forces. (b) Equivalent Set of Forces Applied 
to the Element Nodes. (c) Oppositely Acting Forces Required to 
Maintain Static Equilibrium. (d) Static Force Equilibrium 


The nail shear forces F shown in Figure 6(a) can be thought of as 
externally applied loads and can be replaced by an equivalent set of 
forces which act at element nodes t and b (Figure 6(b)). To resist 
these forces and maintain static equilibrium of the "nail element", an 
equivalent but oppositely acting set of forces must be applied to the 
nodes (Figure 6(c)). In matrix notation, these forces are given as: 


Pt F 
tir} =( pp) = -F [4] 
m (Dt+Dp)*F 


where pt = Horizontal force applied to node t. 
Pb = Horizontal force applied to node b. 
m = Sum of the moments applied to nodes t and b, positive 
counterclockwise. 
Dt = Vertical distance between node t and point Cr. 
Dp = Vertical distance between node b and point Cp. 


Since the force F is a function of A, nodal loads pt, pp and m can be 
related to the displacement of nodes t and b once the relationship 
between the nodal displacements and the slip A is determined. To 
relate the displacements of nodes t and b to A, the X displacements of 
points Cy and Cph are first defined in terms of the displacements of 
nodes t and b respectively. That is, 


u of Ce = üt + D_*O [5] 
u of Ch = up - Dp*9 [6] 


where u = X displacement, positive in the positive X direction. 
ut = X displacement of node t. 
up = X displacement of node b. 
Ə = Rotation of all nodes (including t and b) lying on the 
section line, positive counterclockwise. 


In developing the above equations, it is assumed that: 

l. Both nodes t and b are located on the neutral axes of their 
respective members. 

2. Node t and point Cy both lie on a line which is perpendicular to 
the neutral axis of member T, and node b and point Cy both lie on 
a line which is perpendicular to the neutral axis of member B. 

3. The distribution of strains through the depth of the individual 
layers is linear. Thus, sections of a frame member that are 
originally plane, remain plane. 

4. Rotation 9 is small. 


Since the slip A is equal to the difference between the X displacements 
of points Cy and Cy, the following equation can be written. 


A = Urb > Ub + (Dt+Dp)*9 [7] 
or in matrix notation, 


Ut 
A= [1 -l1 DgtDp] id) where {d} =( up [8] 
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In the above equation, A is positive if at the interface between the 
two members, the top member moves in the positive X direction with 
respect to the bottom member. 


: 1 -1 (Dt+Dp) 
(r) = ke*| -1l 1 -(Dt+Dp) | (d) [9] 
(Dt+tDp) = (D_et+Dp ) (Dt+Dp) 2 


By definition, the product of kę and the 3-by-3 matrix is the stiffness 
matrix for a nail element in which the nail joint exhibits linear 
load-slip behavior. It is important to realize that with this 
stiffness matrix and the source code for a conventional plane-frame 
structural analysis program, a new program can be created for the 
analysis of any multi-layered, horizontally nail-laminated assembly in 
which nail joint behavior is assumed linear. 


Nonlinear Nail-Joint Load-Slip Behavior 


As previously noted, the load-slip relationship of a nail joint is 
usually highly nonlinear. This nonlinear behavior must be taken into 
account when the stiffness of a nail-laminated assembly under 
relatively high loads must be accurately predicted. In equation form, 
the nonlinear load-slip relationship of a nail joint can expressed as: 


F = k(A)*A i [10] 


where k(A) is a secant stiffness that is a function of the slip. When 
equation [7] is substituted into thís load-slip equation and the 
resulting equation for F is substituted into equation [4], the 
following element stiffness matrix is obtained: 


1 -1 (Dt+Dp) 
[k(d)] = k(A)*| -1l l - (Dt+Dp) [11] 
(Dt+Dp) -(Dt+Dp) (Det+Dp) ? 


The notation k(d) is used here to indicate that the stiffness matrix is 
a function of the nodal displacements. This due to the fact that k(A) 
from equation [10] is a function of the nodal displacements. 


The form that k(A) takes is dependent on the model used to describe the 
nail joint load-slip relationship. One model that has been used by a 


number of researchers is the three parameter Foschi equation (Foschi, 
1977) which is given as: 


|F] = (Bo + B1*|4]) [1 - exp(-89*|A]/Bo) ] [11] 


where: 


ll 


|F| = Absolute nail shear force as defined by Foschi equation. 
JA = Absolute slip. 

B2 = Initial stiffness. 

8, = Stiffness at large slips. 

Bo = Intercept of the asymptote, with slope 6). 


The plot in Figure 7 of Foschi’s equation illustrates the relationship 
between parameters Bo, 61, and B9. Note that when Foschi’s equation is 
used, k(A) in equations [10] and [11] becomes |Fe|/|A|. 


Fi a 





A 
FIG. 7.- Load-Slip Characteristics of the Foschi Equation 


Another nonlinear load-slip model which has been used by some 
researchers is the following two parameter function introduced by 
McClain (1975). 


[Fm] = £4*log(1 + B5*|Al) [12] 


In addition to the McClain and Foschi equations, models consisting of 
higher order polynomials are commonly fit to load-slip data. 


When the element stiffness matrix is a function of nodal displacements 
the relationship between the loads applied to an element and the nodal 
displacements is written as: 


[k(d)]{d}) = (r) [13] 
Likewise, the set of structural equilibrium equations has the form: 
[K(D)](D} = {R} [14] 


To solve this nonlinear system of equations a technique can be used 
which involves the removal of all nonlinearities from the element 
stiffness matrices before the structural stiffness matrix is 
formulated. The nonlinearites which are removed from an element 
stiffness matrix show up on the right side of the equation in what is 
referred to as an element load modifying vector. 


The removal of the nonlinearites from the element stiffness matrix is 
accomplished by adding the term [k]{d} to each side of equation [13] 
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and then rearranging the equation. 


{k]{d} = {r} + ([k]-[k(d)])(d} 
{k]{d} = {r} + (w(d)} [15] 


where: 


{w(d)} = element load modifying vector, 
= ({k]-[k(d)]){d}, 
[k(d)] = actual element stiffness matrix which is a function of 
nodal displacements, 
[k] = an "approximate" element stiffness matrix. This matrix 
is not a function of nodal displacements. 


As far as the mathematics are concerned, [k] can be any element stiff- 
ness matrix that is not dependent on nodal displacements; however, [k] 
generally represents the initial stiffness of the element or some other 
linear, first order approximation of the element stiffness. For 
example, when the Foschi equation (equation [11]) is used to describe 
the nail joint load-slip characteristics the actual nail element 
stiffness matrix is given as: 


[k(d)] = |[Fe]/la|*{ -1 1 - (De+Dp) 
(Det+Dp) - (De+Dp) (Dt+Dp) 2 


To obtain an "approximate" element stiffness matrix [k] (one that is 
not a function of nodal displacements), |F¢|/|4| can be replaced by 
parameter Bo - the initial stiffness of the nail joint. 


l -1 (Dt+Dp) 
[k] = £2*| -1l l - (Dt+Dp) 
(Dt+Dh) -(Dt+Dp) (Dt+Dp)? 


By definition, the load modifying vector {w} then becomes: 


(Bo - |Fel/lah*a 
{w} = (Bo - |F Fal tc-a) 
(Bo - |Fel/la])*(Dt+tDp)*A 


where A is given by equation [7] and |Ff| by equation [11]. 


Once assembled, element stiffness matrices and element load modifying 
vectors are loaded into the global stiffness matrix [K] and global load 
modifying vector {W(D)}, respectively. Direct iteration (Cook, 1981) 
is then used to solve the structural equation system. The system of 
equations solved by this method has the form, 
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[K] (D) 441 = (R) + (W(Di)) [16] 


where: 


[K] global stiffness matrix. 
{R} = load vector. 
{W(Di)} = load modifying vector which is a function of the 
displacements from the previous iteration. 
{(D} 441 = vector of unknown displacements. 


The global stiffness matrix [K] is formed and reduced only once thereby 
minimizing computing time. All nonlinearities appear in the load 
modifying vector {W}. This vector is calculated using displacements 
from the previous iteration. The load modifying vector is then added 
to the load vector and the equation is solved for the unknown 
displacements (D};4]. Initial displacements are found through solution 
of the linear equation system [K]{D} = {R}. 


OVERVIEW 


The nodal load-displacement relationships developed in the previous 
section(s) form the foundation for a finite element program entitled 
MLBeam that was developed by the author. MLBeam is not an elaborate 
program. With little work, a program similar to it can be developed by 
modifying the source code for any plane-frame structural analysis 
program which incorporates the direct stiffness method of analysis. 


Programs based on the method of analysis presented here can be used to 
analyze a variety of horizontally nail-laminated assembly designs. 
Basically, there is no limit on the number of layers in the assembly, 
nor do the layers have to be continuous or complete. The wood members 
used in an assembly can be of different height, width and MOE. 

However, they should all be located such that they form a symmetric 
assembly with respect to the X-Y plane. Note that no limitation need 
-be placed on the location, number, or variety of fasteners used in an 
assembly. Likewise, no restrictions need be placed on the type, number 
or location of supports. 


If there is a major difference between the method of analysis herein 
described and those embodied in some of the more sophisticated finite 
element programs, it is that some of the more elaborate methods of 
analysis for horizontally laminated assemblies address continuous 
connections between layers (i.e., elastomeric adhesives). When 
continuous connections are considered, terms in the structural or 
global stiffness matrix become more complex. There are also more 
nonzero off-diagonal terms in the matrix. This is due to the fact that 
nodal displacements in the Y direction are not always uncoupled from 
nodal forces which act in the X direction. The fact that some nodal Y 
displacements will be related to X directed nodal forces is illustrated 
in Figure 8. When both nodes on the same end are given a Y 
displacement (with all other nodal displacement degree of freedoms ` 
fixed) there is a certain amount of interlayer slip (although none 
actually occurs right at the member ends). If adhesive is used which 
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resists this slip, X directed forces will have to be applied to the 
nodes to maintain static equilibrium of the free body. 
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FIG. 8.- Interlayer Slip Caused By A Vertical Translation of the Left 
Element Nodes 







Programs like MLBeam can be used to approximate the stiffness of 
assemblies with elastomeric adhesive. This is accomplished by using a 
series of nail elements and assigning each element a stiffness 
equivalent to that which the adhesive has over an area defined by the 
nail element spacing. 


SUMMARY 


A nonlinear, two-dimensional, finite element program was developed for 
analyzing multilayered, horizontally nail-laminated beams. Nodal 
load-displacement relationships were formulated for a special nail 
element using direct physical arguments. The nail element is used to 
model nail shear joints between layers. It is the only element in 
addition to the conventional plane-frame element that is needed to 
model horizontally nail-laminated assemblies. 


REFERENCES 


Cook, R. D. 1981. Concepts and Applications of Finite Element 
Analysis. 2nd Edition, John Wiley and Sons, Inc. New York, New 
York. 





15 


Foschi, R. O0. 1977. “Analysis of Wood Diaphragms and Trusses. Part 
I: Diaphragms." Canadian Journal of Civil Engineering, 
4(3):345-352, 





Foschi, R. O. 1982. "Structural Analysis of Wood Floor Systems." 
Journal of the Structural Division, ASCE, 108(ST7):1557-1574. 





Goodman, J. R. 1969. “Layered Wood Systems with Interlayer Slip." 
Wood Science, 1(3):148-158. 


Itani, R. J. and Hiremath, G. S. 1980. "Partial Composite Action in 
Diaphragms." Journal of the Engineering Mechanics Division 





ASCE, 106(EM4) : 739-752. 


Kamiya, Fumio 1985. "Nonlinear Bending Behavior of Nailed Layered 
Beams." Wood and Fiber Science, 17(1):117-131. 





Kuenzi, E. W. and Wilkinson, T. L. 1971. Composite Beams-Effect of 
Adhesive or Fastener Rigidity. USDA Forest Service Research 
Paper FPL152. Forest Products Laboratory, Madison, WI. 


McClain, T. E. 1975. Curvilinear Load-Slip Relationships in Laterally 
Loaded Nailed Joints. PhD Thesis. Colorado State University, 
Fort Collins, CO. 


McCutcheon, W. J. 1977. Method for Predicting the Stiffness of Wood- 
Joist Floor Systems with Partial Composite Action. USDA Forest 
Service Research Paper FPL289. Forest Products Laboratory, 
Madison, WI. 


McCutcheon, W. J. 1984. Deflections of Uniformly Loaded Floors: A 
Beam-Spring Analog. USDA Forest Service Research Paper FPL449. 
Forest Products Laboratory, Madison, WI. 


McCutcheon, W. J. 1986. "Stiffness of Framing Members With Partial 
Composite Action" Journal of the Structural Division, ASCE, 
112(ST7) :1623-1637. 


Segerlind, L. J. 1984. Applied Finite Element Analysis. 2nd Ed., 
John Wiley and Sons, Inc. New York, New York. 


Thompson, E. G., Goodman, J. R., and Vanderbilt, M. D. 1975. "Finite 
Element Analysis of Layered Wood Systems. Journal of the 
Structural Division, ASCE, 101(ST12):2659-2672. 








Wang, C. K. 1983. Intermediate Structural Analysis. McGraw-Hill 
Book Company, New York, New York. 


Wheat, D. L., Vanderbilt, M. D., and Goodman, J. R. 1983. "Wood 
Floors With Nonlinear Nail Stiffness." Journal of Structural 
Engineering, ASCE, 109(5):1290-1302. 





Willems, W. and Lucas, W. M. 1978. Structural Analysis for 
Engineers. McGraw-Hill Company Inc., New York, New York. 


